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Abstract 

A cellular automaton model is used to describe the dynamics of the 
catalytic oxidation of CO on a Pt(W0) surface. The cellular automa- 
ton rules account for the structural phase transformations of the Pt 
substrate, the reaction kinetics of the adsorbed phase and diffusion 
of adsorbed species. The model is used to explore the spatial struc- 
ture that underlies the global oscillations observed in some parameter 
regimes. The spatiotemporal dynamics varies significantly within the 
oscillatory regime and depends on the harmonic or relaxational char- 
acter of the global oscillations. Diffusion of adsorbed CO plays an 
important role in the synchronization of the patterns on the substrate 
and this effect is also studied. 

I. INTRODUCTION 

The study of catalytic oxidation processes on Pt surfaces is of considerable in- 
terest since the apparently simple kinetics of this reaction can give rise to different 
types of behavior as a result of competing surface processes, substrate phase trans- 
formations and global coupling effects. Extensive investigations have been carried 
out in the last decade including both experiment and theory [l],[| . Experiments have 
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demonstrated the existence of a variety of non-linear phenonena such as nucleation- 
induced front propagation, spiral wave dynamics, faceting structures and kinetic 
oscillations. 

In this article we re-examine some aspects of the global oscillations observed in 
the oxidation of CO on Pt(100) surfaces. By "global" we mean oscillations in 
concentrations of adsorbed species or Pt substrate states averaged over the entire 
crystal surface or, in experiments, some macroscopic region of the surface. The 
principal features that give rise to such oscillations, which involve coupling between 
surface phase transformations and the oxidation kinetics, were elucidated in earlier 
studies and are as follows: The surface-catalysed CO oxidation occurs by a 

Langmuir-Hinshelwood mechanism ||: 

CO + *$CO ads , (1) 

2 +2* §20 ads , (2) 

CO ads + O ads % C0 2 T +2* , (3) 
CO + O ads % C0 2 T +* , (4) 

where the * indicates an empty adsorption site. The first three steps (|1|-§|) constitute 
the Langmuir-Hinshelwood mechanism. The fourth step (f|) is included to prevent 
poisoning of the Pt surface by adsorbed oxygen, which is not observed in experiment. 
It involves the direct reaction of an adsorbed oxygen atom with CO from the gas 
phase to produce C0 2 in the gas phase plus an empty surface site. 

The actual mechanism of CO oxidation on Pi(100) involves a coupling between 
the reconstruction of the Pt surface and the adsorption and reaction kinetics [|TJ. 
The structure of the topmost layer of clean Pt(100) differs from that in the bulk 
(lxl) and is reconstructed to form a hex structure. || When the CO coverage 
exceeds a critical value this reconstruction is lifted and the surface reverts to the 
lxl form. In general the reaction rates may be different for the two surface 
phases and we have accounted for this possibility by appending the superscript 
a = 1 x 1 or hex to the rate constants in the Langmuir-Hinshelwood mechanism. 
Phase transformations between these two surface forms are triggered by and, in 
turn, influence the characteristics of diffusion and reaction kinetics. It is the switch 
between two surface forms with high and low reactivities that is responsible for the 
existence of global kinetic oscillations. 

Apart from these general observations on the origin of the kinetic oscillations, it 
is known that for Pt( 100) fllU|the spatial coupling of different parts of the surface is 
caused by the propagation of chemical fronts of oxidation across the crystal surface. 

The nature of the wave propagation processes depends on both the type of 
global oscillation and the diffusion coefficient of adsorbed CO, D CO - We investigate 
the spatial structure that underlies the global oscillations. 



2 



The most extensive theoretical descriptions of wave propagation processes associ- 
ated with CO oxidation on various Pt surfaces have been based on reaction-diffusion 
equation descriptions of the spatiotemporal dynamics fLl] , but lattice-gas cellular 



automaton models @J^,|T3 have also been employed in such studies. These lattice- 
gas automata are probabilistic cellular automata that employ a particle description 
of the dynamics and should be distinguished from other cellular automata that are 
designed to simulate reaction-diffusion equation dynamics. |I3|-|T5[| While most of 
the phenomena of interest occur on sufficiently long space and time scales that a 
continuum description using reaction-diffusion equations is appropriate, the cellular 
automaton models have the advantage that they naturally incorporate fluctuations 
that are responsible for the nucleation events that play a role in the process and, as 
well, can be extended to smaller scales where continuum models are no longer appro- 
priate. Cellular automaton models allow one to describe the system at a mesoscopic 
level which incorporates the essential mechanistic features of the reaction dynamics. 
Several such models for surface catalysed CO oxidation employing simple Langmuir- 
Hinshelwood kinetics have been developed in the past years |TB|,|T7[. 

In this paper we use a cellular automaton model for catalytic oxidation of CO on 
Pt(100) which accounts for the coupling among the local reaction kinetics, substrate 
phase transformations and diffusion of adsorbed species. This model is described in 
Sec. y. In Sec. [TTJ the model is used to explore various aspects of the spatiotem- 
poral dynamics in the regime of global oscillations and, in particular, to investigate 
the spatial structure that gives rise to these oscillations. Section [TV] contains the 
conclusions of this study. 



II. CELLULAR AUTOMATON MODEL FOR CO OXIDATION 

From the description given above we see that any model of this surface catalysed 
reaction must account for the coupling among phase transformations, diffusion and 
reaction. The cellular automaton model we employ has features that are like those 



of the Monte Carlo models [If:] in that asynchronous site updating rules are used, 



but there are also features in the model that are adopted from reactive lattice-gas 



automata [18] and, for the surface phase transformation process, rules similar to 



those of traditional cellular automata are used. |T9| The dynamics is described by 
the following rules: 

The cellular automaton dynamics is assumed to take place on a set of regular two- 
dimensional lattices C T . One lattice, C s , is used to represent the Pt(100) substrate 
in its various phases (1 x 1 or hex). The dynamics of the adsorbed CO and O species 
take place on the C a lattice. Each node of C a can support only one particle, either 
O or CO. The nodes of these two lattices have identical labels r. 

Adsorbed species at node r on C a can react with molecules or atoms on neigh- 
boring sites r' G A/" a (r) where A/" a (r) is some neighborhood of the node r. Similarly, 
we let A/" s (r) be the neighborhood of node r on the substrate lattice C s which enters 
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in the description of the surface reconstruction process. Also, we define A/" r (r) to be 
the node r plus its neighbors. In the present work we consider square lattices with 
a total of N nodes and take A/" T (r) to be the Moore neighborhood consisting of the 
eight nearest and next-nearest neighbors of the node r. 

The macroscopic rate constants defined in (|I]-fD are used to construct the micro- 
scopic reaction probabilities. We define an additional rate constant kt r that gauges 
whether a phase transformation will be attempted. For each such reaction rate con- 
stant ki we define the corresponding probability pi = ki/k where k = Thus, 
in the automaton rule described below Pi will specify the probability with which any 
reaction in steps or phase transformation are carried out at a lattice node. 

There are three elements in the cellular automaton rule corresponding to the 
physical processes described above: reaction (to simplify terminology we call ad- 
sorption and desorption "reactions"), phase transformation of the Pt substrate and 
diffusion of adsorbed species. Each automaton time step consists of iV applications 
of the composition of the following operations: 

A node r is picked at random. For node r a reaction i, including phase transfor- 
mation, is chosen randomly with probability pi. 

Reaction 

If one of the reaction steps (|l]-fD is selected and the corresponding z-th reaction 
is sterically allowed, the reaction takes place. In the case of reactions involving 
simultaneously two nodes (for example, O adsorption) the second neighboring lattice 
node is chosen randomly from A/" a (r) and the reaction is carried out if it is sterically 
allowed. The reaction probabilities account for the overcounting of identical pairs 
of nodes. 

Phase transformation process 

If the phase transformation step is selected then substrate phase transformations 



occur by the following mechanism ||12|| . The substrate can exist in two phases, lxl 
and hex. Consider the node r on C s and let n^ (r) = Er'eAT a (r) 8co( r> ) be the 
number of CO molecules at node r plus its neighborhood. Here #co( r ) = 1 if node 
r is occupied by CO and zero otherwise. If we let n be the number of neighbors in 
M a then pg (r) 

— n co( r )/( n + 1) i s the fractional coverage of the region consisting 
of the node r and its neighborhood J\f a . 

• 1 x 1 domain formation mechanism: 

if all nodes in Af s (r) are in the hex phase and A/" a (r) in C a is covered by CO, 
then the node r is changed to 1 x 1 with a probability kjji- 

• 1 x 1 domain growth mechanism: 

if node r on C s is in the hex phase and Pcoi v ) > Pmax then node r is changed 
to 1 x 1 with a probability F{n\ x x)kGi, where n lxl is the number of nodes in 
A/" s (r) in the lxl phase. Here p m ax is some lower fixed fractional coverage 
necessary for 1 x 1 domain growth. 
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• hex domain formation mechanism: 
if if all nodes in J\F s (r) are in the lxl phase and Pco( r ) = then node r on 
C s is changed to hex with a probability km- 

• hex domain growth mechanism: 
if node r on C s is in 1 x 1 phase and if pco( r ) < Pmini then node is r changed 
to hex with a probability F(nh ex )kGH, where rihex is the number of nodes in 
A/" s (r) in the hex phase. Here p min is some fixed upper fractional coverage 
necessary for hex domain growth. 

We have introduced a function F{n a ) {n a > 0) of the neighborhood occupancy, 
= J2 r ' #a( r ')> ( r ' e -Ms( r )) a — hex, 1 x 1), where 6 a (r') = 1 if r' is in phase a 
and zero otherwise. This function models the local structure and is responsible for 
the phase transformation process. It can be used in a phenomenological fashion to 
change the characteristics of this process. The form of F depends on the microscopic 



n 



details of the phase transformation process. | 20| 



Diffusion 

No diffusion in the substrate is permitted. On C a diffusion occurs through random 
jumps of the adsorbed species at node r to one of the empty nodes in J\f a - Some 
flexibility is allowed in the number of diffusion steps compared to reaction steps. One 
may pick additional nodes at random and carry out diffusive jumps. In this way large 
diffusion coefficients can be simulated by increasing the number of diffusion steps 
per reaction step. Similarly very small diffusion can be simulated by increasing the 
number of reaction steps for single diffusion step. A convenient variable to gauge the 
ratio of diffusion to reaction is 5 = (number of diffusion steps) /(number of reaction 
steps). In addition, the diffusion rule can be modified to account for the fact that 
the diffusion rates may be different on the hex or 1 x 1 phases. Thus, the probability 
of a jump from node r may depend on the substrate state at that node. 



III. SPATIAL STRUCTURE AND GLOBAL OSCILLATIONS 

In order to study how reaction, diffusion and phase transformation can lead to 
global oscillatory behavior in CO oxidation on Pt(100) and to examine the un- 
derlying pattern formation that may accompany it, we have carried out a series 
of simulations of the automaton. In the results presented below we have taken 
k\ xl = k\ ex , kif = k hex = 0.001, kf x = 0, k]_f = k h _ e £ = 0, k\ xl = kf x = 0.25, 
k\ xl = k% ex = 0.001 and kt r = 0.6. These values reflect the fact that the major rate 
constant differences between the hex and lxl phases occur in the 2 adsorption 
on the surface: the sticking probability on the hex phase is much smaller than that 
on the lxl phase. We have used k\ xl and k\ xl as bifurcation parameters. In the 
implementation of the phase transformation step we suppose F(n a ) = v? a . The pa- 
rameters that enter in the phase transformation steps are p max = 0.89, p m i n = 0.11, 
k\n = km = 0.001, k GH = k G1 = 0.05. D C o was varied by changing the value 
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of the parameter 5 defined in Sec. |I|. Typically CO diffusion is much larger than 
that of and for simplicity in this study we have taken D^f x = D Q xl = and 
D c £ q = D c x l = D C o and investigated the spatial and temporal structure as a 
function of Deo an d the kinetic parameters. 

Earlier investigations have shown that global oscillations on Pi (100) are 
found in a cusp-shaped domain in the (pco,Po 2 ) parameter plane. In our model 
these partial pressures are incorporated in the kf and k% rate coefficients. We have 
confirmed that our model reproduces this structure in the {k\ xl ,k\ xl ) parameter 
plane and, furthermore, have examined how this bifurcation diagram ETJ] changes 



as Dqo varies. These results are presented in Fig. [I] which shows the region of 
oscillation for two values of the diffusion coefficient as gauged by 5; (a) is a large 
value of Deo with 8 = 10 while (b) corresponds to a small value with S = 1/2. 
The major qualitative feature one observes is that the oscillatory region persists to 
higher k\ xl (pco) values when D C o is small. An understanding of these "phase" 
diagrams requires an examination of the spatial structure that gives rise to the 
global dynamics. We now characterize the global oscillations and then discuss the 
underlying spatial structure. 



A. Global oscillations 

It is difficult to precisely classify the nature of the bifurcations that are responsi- 
ble for oscillatory behavior since our cellular automaton model is a high-dimensional, 
particle-like model and we have no recourse to a simple description based on ordinary 
differential equations (ode). Nevertheless, some general features of the bifurcation 
structure are consistent with those for the global oscillations on Pi(110) described 
in Ref. ]22| where a two-variable ode model that did not include the details of the 



spatial structure was studied. Figure ^ presents phase space plots of the total num- 
ber of COads, O a ds molecules on the surface, Nqo and No, respectively, and the 
total number of sites in the hex phase, N hex . (Clearly, A lxl = N — N hex .) These 
three variables define the phase space (N C o, N Q , N hex ) for spatially averaged dy- 
namical variables. This figure shows the phase space trajectory for these globally 
averaged fields for fixed fc 2 x = 0.04 and varying k\ xl , moving from low to high 
values of this parameter. The results suggest that the oscillations arise by a su- 
percritical Hopf bifurcation on the leftmost boundary of the oscillatory region since 
the (noisy) global limit cycle is harmonic in character and its amplitude grows as 
one moves away from the boundary towards higher k\ xl (cf. panels (a) and (b) in 
this figure). This noisy periodic attractor continues to grow in phase space size and 
evolves to the relaxation oscillation depicted in (c) close to the rightmost bifurcation 
boundary. In this parameter region there are distinctly identifiable slow and fast 
portions to the cycle. Note also that periodic attractor is far less noisy than in the 
harmonic regime. When the bifurcation boundary is crossed this large noisy global 
limit cycle suddenly crashes to a fixed point solution for the global concentration 
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variables. Such a sudden collapse to a fixed point solution was found in ode models 
of oxidation on Pt(llO) |22|. 

Figure ^| shows N hex versus t for a small value of D C o (8 = 1/2) for two values 
of k{ xl along a line k\ = 0.04 in Fig. |l](b). These two values lie near the two 
bifurcation boundaries which showed harmonic and relaxation oscillations for large 
Deo- The global oscillations are highly irregular and much less well resolved, re- 
gardless of the location within the oscillatory regime. We have chosen to plot only 
the time variation of a single variable rather than a phase plane plot since the irreg- 
ular character of the oscillations leads to a complex phase space trajectory where 
the oscillatory structure is obscured. As a result the bifurcation boundaries shown 
in Fig. |I|(b) are subject to large uncertainties. 

B. Spatial structure 

We now consider the spatial structure that underlies these global oscillations. 
While it has been established that chemical waves propagating across the surface 
are responsible for the existence of the oscillations, we shall show that the character 
of these surface structures varies considerably depending on the region of parameter 
space within the oscillatory domain of the bifurcation diagram and the magnitude 
of Deo- 

The most regular global oscillations are found for large Deo in the oscillatory 
region of the bifurcation diagram corresponding to high k\ xl where one observes 
relaxation oscillations. The spatial pattern formation that accompanies these oscil- 
lations is shown in Fig. |] where the states of both the Pt surface and the adsorbed 
phase are shown for a series of six times within one global oscillation period. The 
upper panels in each time group (labeled Pt) show the Pt surface phases while the 
lower set of panels in each time group (labeled ^4) shows the state of the adsorbed 
phase. The Pt surface exhibits well defined propagating fronts that correspond to 
conversion of the Pt surface between the hex and lxl phases. At the begining 
of the cycle shown (upper left) most of the surface is in the hex phase (black) and 
there are islands of the lxl phase (white). Since the CO coverage is high the more 
stable surface state is the lxl phase and the lxl islands quickly grow and consume 
the hex phase until the entire surface is 1 x 1 with a few very small islands barely 
discernible in the top right panel of the figure. With the surface in the lxl phase, 
O adsorption is high, CO is consumed and its coverage falls. The hex islands now 
grow and consume the lxl phase until the surface reverts to a nearly complete 
hex structure with a few small lxl islands. If one observes the spatiotemporal 
dynamics over many global oscillation periods one sees a constantly changing pat- 
tern of waves of phase transformation sweeping over the surface but the dynamics 
is sufficiently well synchronized that the cycling time between the "/iex-rich" and 
"1 x 1-rich" phases is nearly constant. 

In the model the adsorbed O atoms do not diffuse but CO does. Apart from 
some small scale fluctuations the adsorbed species are uniformly distributed and do 
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not exhibit chemical waves (cf. A panels in Fig. [|). In this case Deo is sufficiently 
large (5 = 10) that CO diffusion is able to maintain homogeneity over the finite 
simulation domain. A rough estimate of the diffusion length is £jj ~ (5t osc Y^ 2 , where 
t osc is the oscillation period. For the conditions of the simulation described above 
we have t osc = 640, £r> ~ 80 and N = 100 so that the ratio of the diffusion length 
to the system size is £d/L ~ 0.8 and one expects diffusion to homogenize the CO 
concentration gradients, as is indeed observed. The adsorbed species concentrations 
oscillate but there is almost no pattern formation in the adsorbed layer. 

When Deo is smaller (5 = 1/2) there are time episodes where the waves of phase 
transformation that sweep the surface are no longer synchronized. This desynchro- 
nization is evident in the plots of Nh ex versus t shown in Fig. [5] where one observes 
trajectory segments with large amplitude and others with small amplitude. In the 
time segments where the global concentrations have small amplitude the spatial 
structure is highly disorganized. Figure || shows the spatial structure during such 
a small-amplitude episode where one sees that the irregular time variation of the 
global trajectory arises from the irregular surface structure: while waves of phase 
transformation still sweep across the Pt surface, the nucleation and growth events 
have a highly random character. The surface does not revert to a nearly pure 1 x 1 or 
hex form since new nucleation and growth events occur before others are completed. 
For the system parameters in this case we have t osc ~ 1100 and £o ~ 24 which is 
much smaller than the system length, £d/L ~ 0.24. In this circumstance CO dif- 
fusion will no longer be able to maintain spatial homogeneity in the concentration 
across the surface. This is evident in the lower panels (A) in Fig. [| which show 
the local concentrations of the adsorbed species: there are now significant spatial 
variations in the adsorbed species concentrations. 

Of course, the general structure of these results is to be anticipated on physical 
grounds. A fixed value of the diffusion coefficient and kinetic parameters implies 
a fixed diffusion length £&. As the system size increases it will be impossible for 
diffusion to maintain a synchronized dynamics across the entire surface and one 
should observe phase turbulence. [^3| If the results for small diffusion are viewed 



on smaller length scales, small compared to the diffusion length, then synchronized 
dynamics will again be observed. The domain over which the global concentrations 
are computed in experiment is related to the area of the surface probed by the 
various measuring techniques. This point has been discussed in Refs. 0,0. Apart 
from these rather general observations it is interesting to recall that in the present 
model only CO diffuses and it is the Pt substrate state that exhibits the most well 
defined surface structure. This structure is best organized when there is no spatial 
pattern in the adsorbed phase. In this sense the standard arguments leading to phase 
turbulence take a somewhat different character that makes the results presented here 
less obvious. 

Next, we consider the behavior near the bifurcation boundary where the oscilla- 
tions are harmonic in character and appear to arise from a Hopf bifurcation. Here 
the situation is quite different, even for the large diffusion case where the global 
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oscillations are well defined. In view of the fact that the global oscillations grow in 
amplitude and are small just within the bifurcation boundary, we need to examine 
how it is possible for such small- amplitude oscillations to arise in the spatially dis- 
tributed medium. This can be understood from a study of Fig. ||| which shows the 
state of the Pt surface during one global oscillation period. The open squares on 
the curve showing Nh ex versus t in lower panel of this figure indicate the times in 
the global oscillation cycle corresponding to the Pt surface structure in the upper 
panels. In contrast to the rapidly growing waves of phase transformation that sweep 
across the surface in the relaxation oscillation regime, here one observes that large 
portions of the Pt surface remain in the hex phase throughout the global oscillation 
cycle. The growth of hex now largely occurs by repeated nucleation events which 
tend to "fill in " the lxl seas. The hex islands do, of course, grow and shrink 
to some extent but one does not observe the well-defined waves of phase transfor- 
mation as in the global relaxation oscillation case. One sees that some hex islands 
persist throughout the global oscillation cycle and their shape changes slowly over 
several cycles. Consequently, a different spatial mechanism underlies these harmonic 
oscillations. 

IV. CONCLUSION 

The automaton described in this paper incorporates the reactive processes in 
the adsorbed layer as well as phase transformations in the substrate. We have 
demonstrated that it is able to capture the essential of the features of catalytic CO 
oxidation oscillations on Pt(100). 

The main emphasis of this study was the exploration of the spatial structure that 
underlies the global oscillatory dynamics. The nature of the global oscillations was 
found to depend sensitively on the coupling between the substrate phase transfor- 
mations and the dynamics of the species in the adsorbed layer. This is evident from 
the variations in the oscillatory structure that occur when the diffusion coefficient 
of adsorbed CO is varied. The oscillations are regular in character when Deo is 
large, or irregular, and eventually destroyed, when Deo is small. This behaviour is 
connected with the formation and synchronization of hex and lxl islands on the 
substrate. In the case of large diffusion the hex and lxl islands are large but the 
adsorbed layer is disordered. When CO diffusion is small the islands are smaller 
in size but exist in the substrate and adsorbed layers. The dynamics in these two 
layers are strongly synchronized in this case. 

The form of the spatial structure in the Pt surface also depends on the har- 
monic or relaxational character of the global oscillations, which can be tuned by 
varying pco, the partial pressure of CO. In the relaxation oscillation regime well 
defined waves of phase transformation sweep over the surface. One wave passes 
before new nucleation events begin another cycle the phase transformaion process. 
Consequently one sees very regular global oscillations. In the harmonic region, close 
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to the bifurcation boundary, the global oscillations have small amplitude. As a re- 
sult not all the surface is converted to the 1 x 1 or hex phases. There are no well 
defined waves of phase transformation, instead certain islands remain locked in one 
surface phase over several global cycles and nucleation events fill the "seas" between 
the islands. A rich surface structure underlies the global temporal evolution. 

This automaton is a promising tool for study more complex surface catalytic 
oxidation reactions. Both the adsorbed layer kinetics and the nature of the surface 
reconstruction can be modified from that considered here. In particular, detailed 
studies of pattern formation at the mesoscopic level PUS should be possible us- 



ing such a model and the influence of fluctuations and the possible breakdown of 
deterministic descriptions can be considered. 
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FIGURES 
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FIG. 2. Plots of the system trajectory for 5 = 10 in the (Nco, No, Nh ex ) phase plane 
for k\ xl = 0.04 and (a) k\ xl = 0.005, (b) k\ xl = 0.006 and (c) k\ xl = 0.0167 
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FIG. 4. Evolution of the Pi(100) substrate and adsorbed phase in the regime of regular 
relaxation oscillations for k\ xl = 0.04, k\ xl = 0.015 and S = 10. On the Pt surface labeled 
Pt the black areas correspond to the hex phase and the white areas to the lxl phase. The 
color coding for the species lattice labeled A is as follows: black for O, gray for CO and 
white for empty sites. Time increases from t\ on the upper left to tg on the lower right. 
Starting with time t\ the subsequent times are: ti = t\ + 60, £3 = ti + 90, t± = t\ + 420, 
t 5 =t 1 + 450 and t 6 = h + 630. 
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FIG. 5. Evolution of the Pt(100) substrate and adsorbed phase in the regime of regular 
relaxation oscillations for k\ xl = 0.04, k\ xl = 0.017 and 5 = 1/2. Coding same as Fig. |j. 
Starting with time t\ the subsequent times are: = t% + 80, = t\ + 480, £4 = t\ + 840, 
t 5 = t 1 + 960 and t 6 = h + 1110. 
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FIG. 6. Evolution of the Pt(100) substrate in the regime of regular harmonic oscilla- 
tions for k\ xl = 0.04, k\ xl = 0.005 and 5 = 10. Coding same as Fig. |. The lower panel 
shows Nhex versus t. The square points on the indicate the times at which the spatial 
structure of the Pt surface is displayed in the upper panel. 
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